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Abstract To discuss the role of the gelatinous layer (G- 
layer) on the origins of the physical properties peculiar to 
the tension wood fiber (TW fiber), the deformation process 
of an isolated TW fiber caused by a certain biomechanical 
state change was formulated mathematically. The mechani¬ 
cal model used in the present formulation is a four-layered 
hollow cylinder having the compound middle lamella 
(CML), the outer layer of the secondary wall (SI) and its 
middle layer (S2), and the G-layer (G) as an innermost 
layer. In the formulation, the reinforced matrix mechanism 
was applied to represent the mechanical interaction be¬ 
tween the cellulose microfibril (CMF) as a framework 
bundle and the amorphous substance as a matrix skeleton in 
each layer. The model formulated in the present study is 
thought to be useful to investigate the origins of extensive 
longitudinal drying shrinkage, large tensile growth stress, 
and a high axial elastic modulus, which are rheological 
properties peculiar to the TW. In this article, the detailed 
process of the mathematical formulation is described. In a 
subsequent article, some TW properties from a 70-year-old 
Kohauchiwakaede (Acer sieboldianum Miq.) will be ana¬ 
lyzed using the newly developed model. 

Key words Gelatinous fiber ■ Tension wood ■ Cell wall ■ 
Reaction wood • Growth stress 


Introduction 

Tension wood (TW) consists of abnormal tissue called 
gelatinous fiber (G-fiber) because it contains a gelatinous 
layer (G-layer) as the innermost layer of the secondary wall. 
The TW often shows characteristic behavior that is different 
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from the normal wood (NW). A high-tensile growth stress is 
generated on the surface of the xylem in the TW region that 
often becomes ten times as large as that in the NW region. 1 ' 2 
The longitudinal Young’s modulus of the TW becomes sig¬ 
nificantly higher than that of the NW. 1 ’ 3 Furthermore, the 
axial shrinkage in the TW tends to exceed more than 1% 
during water desorption, 3 while that in the NW becomes 
less than 0.5%. Some authors attribute these behaviors to 
the intrinsic properties of the G-fiber. 1,3 

On the other hand, more than a few researchers believe 
that the G-layer is mechanically too compliant to bear a 
large stress generation. They base their argument on the 
facts that the G-layer is often convoluted in the lumen of the 
transverse section that is sampled from the water-swollen 
block, and it can be easily peeled off the lignified layer in the 
same direction during microtoming. This gives the impres¬ 
sion that it is attached only loosely to the remainder of the 
secondary wall. 4 From those observations, they consider 
that the various characteristics of the TW should be attrib¬ 
uted not to the flexible G-layer but to the relatively thicker 
outer layer of the secondary wall (SI), 5 and/or to the rela¬ 
tively thinner middle layer of the secondary wall (S2) whose 
microfibril angle (MFA) is expected to be more or less 
different from that in the NW fiber. 6,7 

It is impossible to adjudge which possibility holds until 
we succeed in directly measuring the mechanical properties 
of the G-layer and the lignified layer which requires isolat¬ 
ing them from each other. However, it may be quite natural 
to consider that the lignified layer in the G-fiber would be 
essentially same as that in the NW fiber, because there is no 
specific difference in anatomical and chemical aspects be¬ 
tween them as pointed by Okuyama et al. 1 Thus, the author 
expects that characteristic behaviors peculiar to the TW 
may be attributed to the G-layer. 

Simulation using a wood fiber model is one of the most 
effective methods to estimate the internal properties and 
fine structures of each constituent material in wood cell 
walls. 8,9 In our previous studies, we developed a wood fiber 
model consisting of the compound middle lamella (CML), 
SI and S2. A basic formula was derived to simulate the 
deformation process of an isolated wood fiber when a 
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certain biomechanical change occurs, such as lignin deposi¬ 
tion, 1011 water adsorption, 1213 or external load induction. 914 
The basic formula contained several parameters that were 
derived from the composite structure of the cell wall layer. 
When comparing the simulated results with the experimen¬ 
tal one, concrete values need to be assigned to the param¬ 
eters in a rational manner. We consider that those values 
reflect intrinsic information on the internal properties and 
fine structure of cell wall constituents. 

The objective of the present study was to clarify the role 
of the G-layer in the origin of characteristic behaviors of the 
TW, and then, to improve our previous wood fiber model 12,15 
into a four-layered model having the G-layer as the inner¬ 
most layer. In the present article, the detailed process of 
the formulation on the deformation process of the G- 
fiber model which is induced by a certain biomechanical 
state changes is described. In a subsequent article, some 
TW properties of 70-year-old kohauchiwakaede (Acer 
sieboldianum Miq.) will be analyzed by using the newly 
developed model, including a high-tensile growth stress 
generation, a large longitudinal Young’s modulus, and a 
large axial shrinkage due to water desorption. 


Formulation of the deforming wood 

Wood fiber model 

Figure 1 shows a schematic model of the typical G-fiber 
which consists of the CML, the SI, the S2, and the G-layers. 
Each layer can be approximated as a “two-phase structure”, 
specifically the unidirectional reinforcing element of the 
polysaccharide framework and the encrusting matrix (MT) 
of lignin-hemicellulose compound. The former is mainly 
composed of highly crystallized cellulose microfibril 
(CMF), which is oriented in a certain direction to the fiber 
axis in each layer except the CML. This makes the second¬ 
ary wall and the G-layer mechanically anisotropic. On the 
other hand, in cases in which the orientation of the CMF is 



(a) (b) (c) 


Fig. 1. Multilayered structure of the gelatinous fiber (a), its mechanical 
model (b), and the crosscut surface (c). Each fiber consists of the 
compound middle lamella (CML), the outermost layer of the second¬ 
ary wall (Si), its middle layer (S2), and the gelatinous layer (G) 


randomly distributed in the CML, it is considered to be 
mechanically isotropic. 

In the case of softwood secondary wall, oriented polyose, 
such as (acetyl-) glucomannan, is often arranged along and 
around the CMF, forming the polysaccharide framework 
with the highly crystallized CMF. No oriented polyose has 
been detected in the hardwood cell wall. 16-18 Disoriented 
polyose, that is mainly xylan, is blended with lignin, forming 
the isotropic skeleton of the MT substance. 

In the present study, a single G-fiber is simplified into a 
complex hollow cylinder consisting of the CML, the SI, the 
S2, and the G-layers, as shown in Fig. lb. 

Constitutive equations 

Under low magnification, the reinforcing element of the 
polysaccharide framework is spatially dispersed uniformly 
in each cell wall layer to form the framework bundle. Simi¬ 
larly, the lignin-hemicellulose compound is diffused in each 
layer forming the isotropic MT skeleton. Therefore, it is 
considered that both the framework bundle and the MT 
skeleton occupy the same domain in the macroscopic limit. 
Based on such an idea, Barber and Meylan hypothesized 
the mechanical interaction between the polysaccharide 
framework and the lignin-hemicellulose MT in each layer 
as the following conditions: 12,19,20 

Ojj = ofj + cr“, Sy = ej = e” (1) 

where a 1 , o') 1 . and o tj are the stress tensors in the framework 
bundle, MT skeleton, and cell wall layer as a whole, respec¬ 
tively. e\j, e“, and £ iy , on the other hand, are their respective 
strains. 

In the present study, the deformation of the wood fiber 
model is assumed to be symmetric with respect to the cen¬ 
tral axis (see the condition C4). Thus, we can express the 
fact that torsional deformation of an individual fiber is com¬ 
pletely restricted by the force of binding fibers inside the 
wood specimen. Then a cylindrical coordinate system (O- 
£tr, where €, t, and r represent longitudinal, tangential, and 
radial components, respectively) can be applied to the 
present model as shown in Fig. 2a. 

A visible deformation induced at each point in the cell 
wall is expressed by the observable strain tensor compo¬ 
nents (Ejj, where i,j = £, t, r) in the 0-€tr coordinate system. 
The stress component induced in the MT skeleton of each 
layer (oj“) is related to the observable strain component 
caused in the matrix skeleton (e™, where i, j = €, t, r) as the 
following constitutive equation: 

°ij ~ Cijkl\ £ kl a kl) (2) 

C" kl is an elastic constant tensor, and «") is a tensor of 
internal expansion which is caused by a certain biomechani¬ 
cal change, a”, is observed such that in each layer, a™, = 
e m d kh where d kl is Kronecker’s symbol. £ m is a scalar. C™ t:/ and 
a“, are both assumed to be isotropic. The elastic constants of 
the isotropic MT skeleton (nonzero terms) are denoted as 
follows: 
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a b 




Fig. 2a,b. Coordinate system used in the formulation, a O-ftr local 
orthogonal coordinate system, b Flat board element of the CMF frame¬ 
work bundle. O (in a) is an arbitrary point in the cell wall. Direction of 
x-axis (in b) is parallel to the CMF molecular chains in the SI or S2 
layers 


C m _ m _ y'Tn 

tttt ~ '-'tttt — '-'rrrr 


C m _ y'-im _ f -'m 

Utt ~ '-'ttrr — ^tt tt 


C m _ /~<m _ z -'m 

trtr — '-'rtrt ~ '-'tttt 


l(K + 2S), 

C m _ y'rn _ s~<m 

ttrx ~~ rrtt ~ ^rrtt 



i(x-s). 

(3) 


where K = 3X + 2/t, S = 2u. X and u are Lame’s constants. 

Figure 2b shows a small flat-board element of the frame¬ 
work bundle in the S2 layer, provided that the positive 
direction of normal axis (z-axis) is coincident with the radial 
direction (r-axis) of the fiber model. In this model, the CMF 
and other oriented polyose in the S2 layer are assumed to be 
oriented in an S-helix at an angle of 6, and the one in the SI 
layer is assumed to be aligned normally to the fiber axis. 
The relationship between stress (cr 1 *) and strain (e?) compo¬ 
nents induced in the framework bundle of each layer of the 
secondary wall can be written as the following constitutive 
equation in the O-xyz orthogonal coordinate system: 


olt = Cl **(<£ - alj (4) 

where Cf a % d is the elastic constant of the framework bundle 
in each layer in the O-xyz coordinate system. a'J; is an 
internal expansive strain, caused by a certain biomechanical 
change, that is a diagonal tensor whose components are 


/ f*\ 

/ f f f '\ 

l[aab j = 

= (e ,£ ,£ ) 


where e £ and £ f are the internal expansive strains induced in 
the framework bundle in the directions parallel and perpen¬ 
dicular to the cellulose molecular chain, respectively. 

We supposed that the framework bundle is considerably 
compliant in its transverse direction. Therefore, all shear 
moduli, Poisson’s ratios, and the Young’s modulus in the 
transverse direction are small enough to be neglected. This 
means that the CMF including other oriented polysaccha¬ 
ride cannot be kept in a bundle shape without a reinforcing 


agent. The reinforcing agent is considered to be the lignin- 
hemicellulose isotropic matrix. Then, we may consider the 
stiffness components C f a f cd are all nil except (=E). 

By transforming the coordinate system from O-xyz into 
O-ftr system, Eq. 4 is rewritten into a new expression as 
follows: 

Ofy = - a f kl ) (5) 

where 


Cijid — R ia R< h R kc R Id C, 




na 1 '-jb 1 '-kc 1 \d y ^abed, •> 


a ij — Ri a Rjb a ab 


( 6 ) 


R tj is a transformation matrix between both coordinate sys¬ 
tems. The nonzero terms of the stiffness components of the 
framework bundle {C\ jkl ) are expressed as 


Cuu — c E, 
CL = s 4 e, 
al { = e l c 2 + 


Cun — C utf — C S E, Cuu 
C n tt = C em ~ ~ cs E, Cuu z 
e r s 2 , al = e l s 2 + e r c 2 , a € f , 


Cuu c sE, 
= c 2 s 2 E, 

= — (e f — e f )cs 


(7) 


provided that we considered C\ jkl = C) ikl = C] ilk and C[ jkl = 
C f klij , where c = cos cp, s = simp, and E is Young’s modulus of 
the framework bundle in the direction along the cellulose 
molecular chain. E, S , £ m , £ f , p (=r out /r in ), and cp take respec¬ 
tive values in each layer as assumed in Table 1. They are not 
unknown values to be solved but known constants to be 
given in advance. 8,13 S is the shear modulus (X2) of the 
matrix skeleton in each layer. In the same way as E, S takes 
a respective value in each layer. S is denoted as S 0 , S u S 2 , and 
S 3 in the CML, SI, S2, amd G-layers, respectively. 

By rearranging Eqs. 1 using Eqs. 2, 3, 5, and 7, and con¬ 
sidering the compatibility of strains for the axisymmetrical 
deformation without torsion, 


£ r - £, 



we obtain 


a t = -Ks m - Ec 2 e f + + 2 S) + Ec^e ( 

+ l* { K-S) + E C ‘f}', + l[K- S y£, 

o t = -KE m - Es 2 e 1 + ||(A: - 5) + EcV 

+ {i(2X + S) + & <} e , + i(X- S vf, 

°i = —Ke m +|(A: - % +|(2 K + S)E t + ±(K + 2 S)r^, 
a u = — Ecs(e ( c 2 + E t s 2 — £ f ), = o T( = 0 

( 8 ) 


where a, = o, h a t = a tt , a r = a IT . s ( = e u , £ t = £,„ £ r = £ rr . r is 

a radial distance from the central axis. Among the shear 
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Table 1. List of the parameters in the basic formulae A and B 


Layer 

<P 

r in 

^out 

P 


5 

E 

Pin 

Tout 

L 

£ m 

£ f 

CML 

- 

n 

r 0 (=r, + h) 

Pa 

h: 

So 

E 0 (=0) 

Pi 

Pa (=0) 

L 0 

0 

0 

SI 

90 deg 

r 2 

n 

Pl \ 

H) 

St 

E 3 

Pa 

Pt 

Et 

c m 

£ 1 

J 

£ 1 

S2 

e 

r 3 

r 2 

Pa 

H] 

| s 2 

e 2 

P 3 

Pa 

Li 

c m 

£ 2 

J 

£ 2 

G 

Odeg 

U 

r 3 

Pa\ 

H) 

s 3 

e 3 

P*(= 0) 

Pa 

L 3 

P m 

£ 3 

J 

£ 3 


qj, the microfibril angle; r ta , inner radius of each layer; r out , outer radius of each layer; p = r out /r in ; S, shear modulus (X2) of the matrix skeleton; 
E, Young’s modulus of the framework bundle in the direction along the cellulose molecular chain; P to , inner pressure; P out , outer pressure; L, see 
Eqs. 10; e”, inelastic strains in the matrix skeleton; e £ , inelastic strains of the framework bundle in the direction parallel to cellulose molecular 
chains; CML, compound middle lamella; G = EIS ; G, = EJS{, G 2 = E 2 /S 2 ; G 3 = E 3 IS 3 \ Q = FlEp, F = S a -h!r ] ; M = S 2 IS 3 ; N = S 3 /S 2 


stress components, only a (t is not null, which is an inevitable 
consequence from the assumption of the axisymmetrical 
deformation, i.e., s tt = 0. In this study, stress equilibrium for 
axisymmetrical deformation was assumed as follows: 


Deriving the basic equations 
Assumptions 

The dimensional change of the wood fiber model can be 
expressed as a set of normal strains, namely, e e in the longi¬ 
tudinal direction, and £ t (r)| r=v £ t (r)| r=v £,(r)U 2 , £,0)l v 
and £ t (f)| r=r in the tangential directions at the respective 
radius. In the same way as described in previous work, 1112 ' 14 
based on Eqs. 8, we solve e ( , £,(r)U v £,(/)U , £,0)1,. v 
£,(r)| r=r , and £ t (r)| r=r under the conditions Cl, C2, C3, and 
C4, and assumptions Al, A2, and A3: 

Cl. E e is constant for all r. This is based on the assumption 
that the wood fiber model is an infinitely long cylinder in 
the 1-direction. We denote e t as £ L hereafter. 

C2. 0,0)1, , 4 = -P 4 , K(Ol, r n = -P 0 - 0, = O = 0. These 
boundary conditions mean there are no internal and ex¬ 
ternal pressures acting on the wood fiber model (see 
Table 1). 

C3. The external force (/ J ,j induced parallel to the wood 
fiber model satisfies the following equation: 

P L = 2kL = J o,dA = J op'drdO = 2n^ro,(r)dr 

crosscut crosscut 

surface surface 

( 10 ) 

C4. The deformation of the wood fiber model is assumed to 
be symmetric with respect to the central axis. 

In addition to the conditions, the assumptions Al, A2, 
and A3 were made at processing: 

Al. K 2>> 5, so that SIK is enough small to be negligible in 
each layer. This hypothesizes that the bulk modulus of 
the MT skeleton ( K ) is larger than 5. This postulates 


that the Poisson’s ratio of the matrix skeleton is almost 
0.5, similar to a kind of elastomer. 

A2. It is considered that E, S. £ m , £ f , and p (=r out /r in ) tend to 
change their values during a certain biophysical change; 
however, those changes can be neglected in the case 
that biomechanical changes were infinitesimally small. 
A3. E = 0 in Eq. 7 in the CML. Namely, E 0 = 0. This does 
not mean that there is no oriented polysaccharide 
framework in the CML, but means that mechanical 
contribution of the randomly distributed polysaccha¬ 
ride framework in the CML should be isotropic. For 
convenience, we assume that 5„ means the shear modu¬ 
lus (X2) of the CML itself. Moreover, (fi/r,) 2 — 0, £„ = 
0 and £™ = 0 in the CML. These conditions assume that 
h (thickness of the CML) is smaller than the SI and the 
S2 layers. 

Basic formula A 

Solving the second and third formula in Eqs. 8 for £ t and 
rdsjdr , and eliminating the term of o t using the equilibrium 
condition of Eq. 9, we obtain 

[25 + + u(S + 2£x 4 )]e, = 35 ■ f m + (1 + 2 u)E ■ s 2 ■ £ f 

- [5 + Ec 2 s 2 + u(-S + 2EcV)]£ L + (1 + 2 u)r^- 
+ 3uo r , 

[25 + Es 4 + u(S + 2Es 4 )\r— =3 E ■ s 4 ■ e m 
—(2 + u)E ■ s 2 ■ £ f + (2 — 3 s 2 + iijE ■ s 2 • £ l 

— (2 + u)r -^ L + 3 Gus 4 o t 
V ’ dr 

(ID 

where u = SIK and G = EIS. Moreover, combining them by 
eliminating the term of £„ we obtain a differential equation 
of r: 
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Table 2. List of the coefficients in the basic formulae A and B 


Layer 

Location 

r 

A 

X 

A 

Q 

0 

z 

CML 

r = r 0 

H, (=i) 

a; (= 2 ) 

G' (= 3 ) 

G(= 0 ) 

Q 0 ( = l) 

< 2>0 (= 1 ) 

G( = 0 ) 


r = r, 

G(= 1 ) 

A (= 2 ) 

G, (= 3 ) 

4 >(= 0 ) 




SI 

;• = r, 

n 


x; 

g 

q 3 


G 


;• = r 2 

c 



A 




S 2 

r = r 2 

n 

^2 

G 

^2 

q 2 

0 2 

G 


r = r 3 

r 2 

^2 

G 

A 2 




G 

r = r 3 

G (=i) 

2 l' ( = 2 ) 

XI (= 3 ) 

G(=o) 

q 3 

0 3 ( = 1 ) 

G 


r = r t 

G (= 1 ) 

A 3 (= 2 ) 

X 3 (= 3 ) 

a 3 (= 0 ) 





(1 + 2»)r 2 d °f + 3(1 + 2u)r^ 1 - — 3 Gus 4 o t 
= 3E ■ s 4 e m - (2 + u)E • s¥ + E • s 2 (2 - 3s 2 + u)t: { 
General solution of Eq. 12 is given as 

1 


( 12 ) 


o T = cy 


C 7 r 


—[3£ 


1 + 2 u 2a L 
— (2 + u)E ■ s 2 e f + E ■ s 2 (2 — 3s 2 + ujs L j 


(13) 


X 


1 

( r^ 1 - lY 

1 + a 

l « - 1 J 


where C\ and C 2 are integral constants, and a is described as 
follows: 


a = 1 + 3G ■ s 


1 + 2 u 


= Jl + 3G • s 


K+ 2 S 


Based on the assumption Al, we impose the following 
conditions: 

« —> + 0 then a —> 1 + 0 
which yields 
r al -i 


lim 


lnr 


“-> 1+0 a —1 

Then, we can simplify the solution of Eq. 13 as follows: 

a r = C,' + C 2 r 2 + | [3 E-s 4 E m - 2E-sV 
+ (2 - 3s 2 )fi • s 2 £ L jlnr 

Under the boundary conditions \o r {r)\ r r = - / J m and 
\a,(r)] r r = ~P olt , we can decide the integral constants C[ 

out 

and C 2 , and we obtain the following solution: 


o r = ~P:„ - 


P 2 - rljr 2 


Pout - Pin + 


4„m 
S £ 


p 2 - 1 f oul ~ ra 2 1 

- 2 E • sV + E ■ s 2 (2 - 3s 2 )e L ]lnp| + |[3 E ■ sV 

- 2E ■ s 2 £ i + E ■ s 2 (2 - 3s z )e L ]ln- 


r 

r in 


( 14 ) 


Then, we substitute the above solution into the first part of 
Eq. 11, and we obtain basic formula A: 


r ■ S L + A- £ t = X ■ £ m + A- £ l 


1 

(rL] 

Pin 

( pI ~ i) 

l r2 J 

S 


1 


p 

1 out 

U - !) 

U 2 J 

5 


where 


r = 1 + i-G-s 4 + G • s 2 (2 - 3s 2 ) 

3 


lnp 


( 2 1 


U - 1 ) 


Xj 


A = 2 + Gs\ X = 3 + -G s 4 - 3G s 4 lnp 


(O' - 1) 



( 

r 2 

\ 



'out 


\ 


? 


) 

V 

r 

/ 


and A = 2 G • s 


lnp 

(r 2 \ 

'out 

c ~ i) 

U 2 J 


Coefficients r. A, X , and A are dependent on G, cp, r, r imt , 
and p, provided that these variables take their respective 
values in each layer as shown in Table 2. 


Basic formula B 


We solve Eqs. 8 for £ L and integrate it over the crosscut area 
of each layer. When integrating it, we assume that £ m and £ f 
are independent of r in the respective layers, and £ ( takes a 
constant value (e l ) over the crosscut surface of the cell wall 
because of condition Cl. As a result, we obtain basic for¬ 
mula B: 


Q ■ £ l = O ■ £ m + E ■ £ l + 


2 L A 
3S-r 2 [p 2 -\) 


2 p Eput I y^1 _ 2^_ 1 

3 (p 2 - l) S L 3 (p 2 - 


W—f r l 

l) S L lr ='“ 


where 


Q= 1 + -G(l - 3s 2 + 3s 4 ) 


(2 - 3s 2 ) 2 

—T—\np) - i 


vp - 1 J 4 J 
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equations into simpler ones in which the unknown variables 

are e L , e t \ r=v £ t | r= , 2 , £ t U, 3 , and s\ r=r ^. 


0=1- G-s 2 [c 2 - s 2 ) 
+ G 2 • /(2 - 3s 2 ) 


P-1 


ln P - T 


2 = -G(2 - 3s 2 ) - |g 2 -s 4 (2 

L = — J c^rdrdP 

each layer 



1 

4 


Coefficients 12, &, and 2 are dependent on G, <p, and p, 
provided that these coefficients take their respective values 
in each layer as shown in Table 2. 


Formulae to describe the deforming G-fiber 

In the G-fiber model as the whole, twelve equations are 
obtained based on the basic formulae A and B. The equa¬ 
tions are as follows: 


Basic formula A in CML at r = r 0 ( a 0), 

Basic formula A in CML at r = r l (at)'), 

Basic formula B in CML (b0), 

Basic formula A in SI at r = r l (al), 

Basic formula A in SI at r = r 2 (al'), 

Basic formula B in SI (bl), 

Basic formula A in S2 at r = r 2 ( a 2), 

Basic formula A in S2 at r = r 3 (a2'), 

Basic formula B in S2 (b2), 

Basic formula A in G at r = r 3 ( a 3), 

Basic formula A in G at r = r 4 (a3'), 

Basic formula B in G (b3). 


These equations constitute simultaneous equations whose 
unknown variables are 

£ L> £ t|r=r„J £ t|r=iii £ t|r=r 2 ) £ t|r=r 3 ! £ t | r=r 4 ■> 

Po, Pu Pi, P31 P*\ L 0 , Ljy, L 2 , 

According to the assumption A3, Eqs. aO and a0' are degen¬ 
erated by each other, which yields 


Therefore, the unknown variables become 

£ L? £ t|r=o’ £ t \ r-r 2 , £ t|r=r 3 5 £ t|r=r 4 5 

Po, Pl, Pi, Pi, P*\ P«, Pi, Pi, P3 

To solve Eqs. a0'-a3', and Eqs. bl-b3 for the unknown 
variables, the following conditions based on conditions C2 
and C3 are imposed: 


The first equation 

From Eqs. aO', bO, al', bl, a2', b2, b3, cl, c2, and c3, we 
eliminate P 0 , P u P 2 , P 3 , P 4 , L 0 , P 3 , L 2 , L 3 , and L. Thus, we 
obtain the first equation 

fl ll £ L fl 12 £ t|r=/i 3“ fl 13 £ t|r=/i 3“ ^I4£ t |r=r 3 3“ tfl5 £ t|r=r 4 

= b n E -f + b l2 e 2 + b l2> e 3 + Cu4 + c 12 s f 2 + c 13 £ 3 + d n P L 

(15) 

where coefficients a n , a 12 , a 13 , a 14 , a 15 , b n , b 12 , b 13 , c u , c 12 , c 13 , 
and d n are respective functions whose concrete forms are 
composed of 6, p 1? p 2 , p 3 , Q(=FIE l ), G,, G 2 , G 3 , M(=5 2 /5 1 ), 
and N(=S 3 /S 2 ). Detailed shapes of those coefficients are 
described in the Appendix. P L stands for the external force 
induced parallel to the wood fiber model, which is related to 
L 0 , Lj, L 2 , and L 3 as follows: 

P L = 2kL = 2k(L q + Lj + L 2 + L 3 ) = J o e rdrd6 

crosscut 

surface 

= 2 nj°o e rdr 


The second equation 

e L , £,( r] , £ t | r=ri , e t | r v and £,\ r=r4 are unknown variables to be 
solved as solutions of an algebraic equation (Eq. 15). To 
solve Eq. 15 for £ L , £,| r v £,( v £ t | r=r3 , and e\ r _ v there must 
be at least four equations that are constituted by the same 
unknown variables. These equations can be derived from 
basic formula A. 

From Eqs. aO', al', a2', a3, cl, and c2, we eliminate P 0 , P u 
P 2 , P 3 , and P 4 . Thus, we obtain 

a 21 E L 3“ fl 22 £ t|r=/i 3" fl 2 3 £ t|r=/i 3“ U 2 4£ t | r=r3 + 025 £ t|r=r 4 

= b 21 E 3 + b 22 E 2 + b 33 E 3 + C 21 £j + C 22 £ 2 + c 23 e 3 + d 21 P L 

(16) 


The third equation 

From Eqs. aO', al', a2, a3, cl, and c2, we eliminate P 0 , P u P 2 , 
P 3 , and P 4 . Thus, we obtain 

fl 31 £ L 3“ ^32 £ t | r —4 3“ fl 33 £,| r=4! + fl 34 £ t |,.=r 3 + 035 £ t|r=r 4 

— b 31 E 3 + b 32 E 2 + b 33 E 3 + C 3 j£j + C 32 £ 2 + C 33 £ 3 + d 31 P L 

(17) 


The fourth equation 


P 0 = o (cl), P 4 = 0 (c2), L = Lg + L 1 +L 2 + L 3 (c3) 

The unknown variables explicitly required in our study are 
£ l, £ tlr= ri , £ tlf=r 2 ’ £ tlr=r 3 , and £,L, 4 . Thus, by eliminating the 
unknown variables P„, P,, P 2 , P 3 , P 4 , L 0 , L u L 2 , L 3 , and L in 
the following manner, we degenerated the simultaneous 


From Eqs. aO', al, a2, a3, cl, and c2, we eliminate P 0 , P l3 P 2 , 
P 3 , and P 4 . Thus, we obtain 

fl 41 £ L 3“ U 4 2 £ t|r=n 3" fl 43 £,| r= ^ + fl 44 £ t | r=rj + fl 45 £ t | r=r4 

— b 4 1 £| + b 42 E 2 + b 43 £ 3 + £ 4 |£| + C 42 £ 2 + C 43 £ 3 + d 41 P L 

( 18 ) 
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The fifth equation 


From Eqs. aO', al\ a2', a3', cl, and c2, we eliminate P 0 , P u 
P 2 , P 3 , and P A . Thus, we obtain 

fl 51 £ L T U52 e t|r=r! F fl 53 £ t|r=^ F fl54f t |r=/5 F fl 55 £ t | r=r4 

= h 51 £! m + h 52 £ 2 m + h 53 £ 3 m + C 51 £j + C 52 £ 2 + C 53 £ 3 + d 51 P L 

(19) 


where coefficients a 21 -a 55 , b 21 -b 5 


-c 53 , and d t , — c/ 5 , are re¬ 


spective functions whose concrete forms are described in 
the Appendix. 

Equations 15-19 constitute the simultaneous algebraic 
equations whose unknown variables are £ L , £?-(=£j, r ), 
£ r(= £ iU 2 ), < 3 (=eJr=r 3 ). and £ ( 4 (= £ ( U 4 )- The values of e™, s 2 , 
e”, £'[, £ 2 , £ 3 , and P, should be given in advance. From the 
simultaneous equations, the solutions should be expressed 
in the following forms: 


£ L =/ll(p) £ i m +/l2(p) £ 2 m +/i 3 (pK +/i4(pH +/15(P) £ 2 
+ /l6(p) £ 3 + /l 7 (p)^L 

£ t =/2l(p) £ i m +/ 2 2(p) £ 2 m +/2 3 (p) £ 3 m +/ 2 4(p) £ l f +/2 5 (p) £ 2 

+ /m(pK + / 27 (P)P[. 

£ ? =/3l(p) £ i m +/3 2 (p) £ 2 m +/33(p) £ 3 m + / 34 (pK +/^(p)^ 
+ / 3 6(P) £ 3 + / 37 (p)^L 

=/4l(p) £ i m +/4 2 (p) £ 2 m +/43(p) £ 3 m + ^44 (p) £ l + / 4 5 (p) £ 2 
+ /46(P) £ 3 £ + / 47 (P)^L 


d £ l_ = fu (p)^ £ i m + /l2 (p)^ £ 2 m + /l3 (p)^ £ 3 m + /l4 (p)^ £ l f 
+ /i 5 (p)d£ 2 + /i 6 (p)d£ 3 + / l7 (p)dP L 

The Eqs. 20' are divided by dt, and are converted into the 
following differential equations: 

£ L =/u(p)eT +/l 2 (p) £ 2 m +/l3(p) £ 3 m +/l4(p) £ l f +/l5(p) £ 2 

+ fi6 (p) £ 3 + fn{p)PL 

£ t =/2l(p) £ i m +/22(p) £ 2 m +/23(p) £ 3 m + /2 4 (p) £ l £ +/2 5 (p) £ 2 
F f 2 (, (p) £ 3 + 

e? =/3l(p) £ i m +/ 3 2(p) £ 2 m +/33(p) £ 3 m +/ 34 (p) £ 1 £ +/ 3 5(p) £ 2 

+ / 3 6(P) £ 3 £ + / 37 (p)^L 

£ ? =/4l(p) £ i m +/4 2 (p) £ 2 m +/43(p) £ 3 m +/ 44 (pK f +/4 5 (P) £ 2 

+ /46(P) £ 3 + / 47 (P)4 

^ =/5l(p) £ i m +/52(p) £ 2 m + fs 3 {p) £ 3 +/5 4 (P) £ 1 £ +/ 5 5(p) £ 2 
+ /56(P) £ 3 + fs 7 {?)h 

( 21 ) 

where the raised dot notation represents the derivative by t. 
Furthermore, p,, p 2 , and p 3 among the components of p 
should be rewritten into the differential quantities dp t , dp 2 , 
and dp 3 , because they depend on the elapsed time during a 
certain biomechanical change. Due to the assumption of 
axisymmetrical deformation, the following relations are in¬ 
quired among respective layers: 


£ t =/ 51 (pK + fsi{p) £ T + Mp) £ 3 +/5 4 (p) £ 1 £ +/5s(p) £ 2 
F fit, (p) £ 3 + / 57 (p)^L 

( 20 ) 

where coefficients/ n -/ 57 are functions of p, and p is a param¬ 
eter vector whose components are 0, p,, p 2 , p 3 , Q, G l5 G 2 , G 3 , 
M, and N. Among the solutions, £ L represents the strain of 
the G-fiber model in the axial direction. According to 
the condition C4, the wood fiber model deforms 
axisymmetrically so that £?(=£?), e[ 2 , Pfi and Pf are equiva¬ 
lent strains to the diametral deformations at their respective 
radii. These strains are induced by a certain biomechanical 
change in the G-fiber. The thickness of the CML is small 
enough to be negligible if compared to those of the SI and 
S2 layers. Therefore, £[' can be regarded as the strain of the 
diameter in the wood fiber model (=£ T ). 


Pi + dp 3 - 


Pi + dp 2 - 


p 3 + dp 3 = 


r 3 + dr 3 
r 2 + dr 2 
r 2 + dr 2 
r 3 + dr 3 
r 3 + dr 3 
r A + dr 4 


1 + d£ t ri 

1 + dp? 


Pi, 


1 + rfff 
1 + de[ l 
1 + de[ 3 
1 + de[' 


Pi, 


Pi 


( 22 ) 


Because dp 3 de [ 4 , dp 2 ds [ 3 , and dp { dP 2 are regarded as higher 
order infinitesimal quantities, from Eqs. 22, we can obtain 
the following differential forms: 

dp 3 = (rfe? — d£j )p l5 dp 2 = (de? — de? )p 2 , 

dp 3 = (ds? - dtfi )p 3 ( ' 22 ’ 

These equations can be rewritten as differential equations 
of t as follows: 


Developing Eqs. 20 into differential equations 

According to assumption A2, Eqs. 20 are valid under the 
condition that the biomechanical change is infinitesimally 
small enough to be neglected. This means that £ L , e' 1 , P 2 , pf 
and £, 4 in addition to s 3 , t 2 , e 3 , e\, e{, £ 3 , and !\ in Eqs. 20 
should be replaced as differential quantities, i.e., de L , dP x \ 
ds\ 2 , de[ 2 , dPf, dsf, de 2 , dpf, ds[, de\, dej, and dP L . Then, each 
expression in Eqs. 20 should be rewritten as a simple differ¬ 
ential form, e.g., 


p, = (fi? - £[ 2 )p„ i> 2 = (s? ~ s?)p 2 , p 3 = (fi[ 3 - £[ 4 )p 3 

(23) 

Equations 21 and 23 constitute a system of simultaneous 
differential equations whose unknown functions are £ L , eJ 1 , 
£?, pf, pf, p 3 , p 2 , and p 3 . On the other hand, e“, e“, £™, £■(, £ 2 , 
e(, and P L are the functions whose f-dependent shapes 
should be given in advance. 
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Solving the simultaneous differential equations (Eqs. 21 
and 23) 

Solutions 


Equations 21 and 23 give the deformation of the wood fiber 
model, which is induced by a certain biomechanical change 
occurring in the cell wall. In our case, examples of biome¬ 
chanical change are the growth strains (maturation strains) 
during cell wall lignification, which is measured as the re¬ 
leased strains of the growth stress; the swelling strains due 
to moisture adsorption; and the elastic deformation caused 
by external load. 

We assume that change in the physical state in the cell 
wall starts at t = 0 and ends at t = Z. By integrating Eqs. 21 
and 23 from t = 0 to t = Z, we can solve them for £ L , e[', t) 2 , 
£, 3 , £\ 4 , p h p 2 , and p 3 . As the initial conditions at t = 0, we 
assume 

£ L (0) = £[‘(0) = £?(0) = £j 3 (0) = £*(0) = 0, 

Pi(0) = pi, P?(0) = p 2 , P 3 (0) = p 3 

Then, we divide the integral interval into n small parts, and 
denote the integrals of Eqs. 21 in the z-th small interval 
[(/ — 1 )Zln < t < iZIn, i = 1,2,..., n] as A'e L , A‘e[\ 
A‘e r t 2 , A‘e r t 3 , and A‘s r t 4 . For example, A‘e l is calculated as 


i-Zjn 


A', L , J 


ds i 


(i-l)Z/n 

i-Z/n 


-A^dt 

dt 


= J (/n(pK + /12 (p)f 2 m + /is(pR + fu (p# (25) 

(i-l)Z/„ 

+ f \5 (p)( 2 +/l6(p)£s + /l 7 (p)P L )<* 


In each small integral interval, the functions £“(?), £™(f), 
£™(t), e\(J), £ 2 (f), £ 3 (/), and respective components in p, in¬ 
cluding p 1; p 2 , and p 3 , must be given in advance. However, 
p l5 p 2 , and p 3 are unknown functions to be solved from 
simultaneous differential Eqs. 21 and 23. The values of p,, 
p 2 , and p 3 , which should be used in the z-th small interval, 
can be estimated in the following process. At first, we inte¬ 
grate Eqs. 23 in the (z — l)-th small interval. Then, we 
obtain 


(i-l)Z/n 

J 


(i—2)Z/« 




f(i-l)Z/n 

1 

[(i-2)Z/n 




dt 



If n is taken as a large enough number, this equation gives 

/TV = (zrv -zryj.pr 2 (26) 

where p) 2 is equal to p 3 (f) at t = (z — 2 )Zln. Thus, the 
value of p\ \ which is p,(z) in the z-th small interval, is given 
as 

pi 1 = pr 2 + a^p, = (i + zt 1 # - a‘-'f ) ■ pr 2 (27) 

By using simultaneous recurrence equations (Eqs. 25 and 
27), we can integrate simultaneous differential Eqs. 21 and 
23 numerically. 


Thus, the natural strain of the deforming wood fiber 
model at a certain time t (= j ■ Zln, j is positive integer 
smaller than n ) can be derived as follows: 


i 

£ L (r)= lim^zl'f L 

l =1 

j 

t T i*' _ |; ™ 


(t) = limY A l e? 

v ' n->°° 


i =1 



(28) 


Strains £ L (t) and give natural strains of the dimensional 
changes induced in the wood fiber model. However, the 
released strain of the growth stress and the swelling due to 
water sorption are measured as the nominal strains in the 
deforming wood specimen. Then, we need to solve Eqs. 21 
and 23 to give anisotropic dimensional changes of the wood 
fiber model as nominal strains. 

The nominal strain of the deforming wood fiber model at 
a certain time t (= j ■ Z!n) in respective directions can be 
derived as follows: 


a(t) = lim 

' rt—>°o 


ri(i+ 



/3(t) = lim 

' n—> oo 


ri(i+ j#) 



(29) 


If the values of a(Z) and /?(Z) are smaller than 1%, changes 
of p 3 , p 2 , and p 3 are small enough to be neglected. In that 
case, we can obtain the a(t) and /3(t) from the first and the 
second formulae of Eqs. 21, which simplifies Eqs. 29 as 
follows: 

a(?) = limYzl'£ L = f— /3(t) = limYzl'f^ 1 = f— -dt 
n ^~7~i i dt J 0 dt 

(30) 


Examples of simulations 

Generation of the maturation strains due to cell 
wall lignification 

Soon after the deposition of the CMF and the highly ori¬ 
ented polysaccharide framework in the secondary wall of 
the tracheid or libriform wood fiber, the matrix substance of 
hemicellulose and lignin deposits among the gaps of the 
framework bundle. In this process, the wood fiber tends to 
swell or shrink anisotropically, which generates anisotropic 
growth stress in the newly formed xylem. This is because 
free deformation of the individual fiber is restricted to in¬ 
side the actual xylem. This process was formulated by using 
a wood fiber model of the CML, SI, S2, and G layers. In this 
case, irreversible deposition of the MT substance and the 
maturation of the CMF framework are postulated to cause 
a change in the biomechanical state of the tracheid or libri¬ 
form fiber wall. 

In the present study, the free dimensional change of the 
G-fiber is simulated during the wall maturation. i\ should 
be null in Eqs. 21. The integral interval is from the deposi- 
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Table 3. Assumed values for the chemical compositions in each layer of the green G-fiber wall 


Layer 

Cellulose crystal 

Oriented polyose 
(non-crystalline) in the 
poly-saccharide framework 

Isotropic matrix 

Crystallinity in the 

polysaccharide 

framework 

CML 

15 (%) 

0 (%) 

85 (%) 

100 (%) 

SI 

15 

5 

80 

75 

S2 

30 

10 

60 

75 

G 

Case a 

40 

40 

20 

50 

Case b 

60 

20 

20 

75 

Case c 

80 

0 

20 

100 


tions of the CMF framework (t = 0) to the completion of 
the G-fiber wall (t = Z). In this interval, the f-dependent 
patterns of ef, £™, e™, s{, s\, £ 3 , and every component in p 
should be given in advance. As the initial conditions at t — 
0, we adopt Eq. 24. According to the observations, the 
released strains of the growth stress are less than 0.5% both 
in the 7 and the t directions. Then, the simplified formulae 
(Eqs. 30) can be used instead of Eqs. 29 when simulating the 
maturation strains a(f)|, =z and /■ In such a case, p,(/), 
p 2 (f), and p 3 (f) are regarded as approximately constant. 

As the preliminary simulation, we calculated the values 
of a(t)\ t=z (=£ l ) and z (=£ x ) under the conditions of 
the parameters assumed in Tables 3 and 4. The values in 
Tables 3 and 4, except those related to the G-layer, were 
used in our previous simulations on the generation of 
growth strains in the latewood tracheid of sugi ( Cryptom- 
eria japonica). s In the present simulation, the various cases 
for the percentage of the cellulose crystal in the G-layer are 
displayed in Table 3, and the relationship between £ L , £ x 
and was calculated for each case of framework 

crystallinity. 

The result is shown in Fig. 3. According to the observa¬ 
tion in the normal fiber of the softwood or the hardwood 
xylem, the released strain of the surface growth stress be¬ 
comes — 0 . 02 % to —0.06% in the longitudinal direction, and 
0.05%-0.1% in the tangential direction. In the TW xylem, 
on the other hand, the longitudinal released strain often 
becomes several times as large as in the NW xylem, and in 
some species forming considerably thick G-layers, such as 
Robinia pseudoacacia , it exceeds —0.5 %. 2 The present 
simulation explains these phenomena well when the G- 
layer is assumed to shrink considerably in its axial direction 
during the G-fiber maturation. 


Longitudinal elastic constants 


At present, the longitudinal Young’s modulus (E L ) of the 
wood fiber model is calculated based on Eqs. 21 and 23 
under the assumption of the moisture steady state. There¬ 
fore, it is assumed that every component in p is constant. 
Moreover, dd", dc"\ dc"', de[, de l 2 , and de\ are all nil. Then, 
from Eqs. 21, we obtain the longitudinal Young’s modulus 
of the G-fiber (U L ) as follows: 




J_dii 

TTTo dc L 


J_ _ 1_ 

/iv(p) 


( 31 ) 



Fig. 3. An example of the simulation on the relations between 
the growth strains and the values of at t = Z (= T 3 ) in the G-fiber. 
Values of the parameters assumed here are displayed in Table 4. Con¬ 
ditions for the crystallinity of the polysaccharide framework in the G- 
layer (see Table 3), solid lines , case a: dashed lines, case b; shaded tines, 
case c 


The substantial Young’s modulus of the G-fiber can be 
calculated as 

_ 1 dl] _ 1 1 

7r(/o 2 - £ 4 ) de L - l~4 ) /l?(p) 

In this case, an increase in the strain energy induced by 
external load is regarded as the biomechanical state change. 

As a preliminary simulation, the values of E L at t = 7) 
were calculated under the conditions of the parameters 
given in Tables 3 and 4. The result is shown in Fig. 4. In the 
case in which the G-layer contains a certain amount of 
crystalline cellulose, calculated values of E L become larger 
with the thickness of the G-layer; however, it is almost 
constant when the G-layer contains no crystalline cellulose. 


Swelling and shrinkage strains due to the moisture 
adsorption 

In the same way as the growth stress case, free dimensional 
change of the G-fiber due to moisture adsorption were 
simulated. Thus, P L should be null in Eqs. 21. The integral 
interval is from the oven-dried state (u = 0 ) to the liber 










Table 4. Assumed values of the parameters for the simulation of the growth strains of the G-fiber 
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Fig. 4. An example of the simulation on the relations between the 
longitudinal Young’s modulus and the value of p 3 in the green G-fiber. 
Conditions for the crystallinity of the polysaccharide framework in the 
G-layer (see Table 3), solid line, case a; dashed line, case b; shaded line, 
case c 


saturation point (u = Z). Therefore, it is regarded that the 
moisture content change u is equivalent to the elapsed time 
t. In the integral interval from u = 0 to u = Z, independent 
patterns of e™, e™, e[, e‘ 2 , £ 3 , and every component in p 

must be given in advance. In this case, the moisture sorption 
in the cell wall causes the biomechanical state change in the 
G-fiber. 

The shrinkage strain of wood cannot be regarded as a 
reciprocal of the swelling strain because two different refer¬ 
ence bases are used in their measurements. If the hysteresis 
effect between shrinkage and swelling processes is small 
enough to be neglected, we can convert the swellings a(ii ) 
and /?(«) into shrinkages a'(u) and /T(i<) at a certain mois¬ 
ture content (u) by using the following formulae: 


a(Z) — a(u) 
a(Z) + 1 ’ 


= p(Z) - flu ) 


P[Z) + 1 


(32) 


In the case that the oven-dried shrinkage of the wood be¬ 
comes more than several percent in the transverse direc¬ 
tion, we should use Eqs. 29 when calculating a(u) and fl(u) 
[or a'(u) and / 3'(u )]. 

In the same way as in the case of growth stress genera¬ 
tion, we can reasonably expect that the longitudinal drying 
shrinkage a' L (u)\ u=z becomes considerably larger with the 
thickness of the G-layer and with the negative value of e 3 £ |„ =z . 


Concluding remarks 

In the subsequent report, further concrete simulations will 
be demonstrated, and the results will be compared with the 
experimental results obtained from the TW of a 70-year-old 
Kohauchiwakaede (Acer sieboldianum Miq.), which are 
already reported . 3 Examples include high tensile growth 
stress generation, large longitudinal Young’s modulus, and 
large axial shrinkage due to water desorption in the G-fiber. 
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Appendix 

Detailed expressions of coefficients fl n -a 55 ; b n -b 53 , c u -c 53 ; d n -d 51 in Eqs. 15-19. 


b u = H- 


aXi_ 

6 M 


{ 1 ) 

fpi-n 

IP? -ij 

1 pi I 


O-A / 

6M 


{ 1 1 

/ 2 1 \ 

Pi ~ 1 

\p2-1J 

2 

A J 


1 V 1(1 

, b u = a d> 2 — -~ta{A ~ 2A), fo 13 = -NaA 2 --- 

6 p 2 2 - 1 

, c 12 = a ■ E 2 — — —a(A 2 — 2/\)„ c 13 = — iVaA 2 | — 2 — 

6 P 2 3 V Pi ~ 1 


y 3 -1 

V P3 


Ps 


J - 2 A 2 
CL ii C( 

3 M 


ft I 1 1 _ ^ ^ , 1 

—- 7 L 2 n ’ “ 21 “ Gl ^ + o' 

p 2 - 1 J 5> 2 2 m: 2 


+ M 

pr J 2 


' p* -i 

P 2 2 


jr 2 + ^MN(p\ - 1 ), « 22 = 2 G,Q, 


f ^2 


pi - 1 ) A _ m 

PI P ’ fl24 2 


A = Imn{p\ - 1), c 21 = i 


^PL_iL 2 + Miv(p 2 _ 1 ), fl25 = 0 , b 2l = ± 
P2 J ^ 


V 

/ _2 


/ _2 


- f 


' p 2 -1 ^ 
P 2 j 


Xi, 


Ptzl U c =" 
Pi 2 J 22 2 


' ft "I 

P2 


I 2 „2 


«3i - GlG + — 


Pi “I 
Pi" 


]a + y(ft - l)A' + |m 7 V ( P 2 3 - 1 ), 


1’ • C-, 0, d 2l 0, 

fl 32 = 2Gj<2, 


( _2 


Pi ~ 1 
pr 


j/f, + y(ft - l)A, «34 = MN[p\ - l), n 35 = 0, b 31 = ^32 = y(ft - l)A, 


fc 33 = |iWiV(p 2 - 1), c 31 = fc 32 = y(p 2 2 - 1)^, c 33 = 0, d 31 = 0, 


«41 = GiG + ±(ft - I)r; + y(p 2 - i)r 2 ' + ^Miv(p 2 - 1), « 42 = 2G 3 g + |(ft - i)/i„ « 43 = y(p 2 2 - i)a 2 , 
fl 44 = MN(pl - l), fl 45 = 0, b 4l = |(p 2 1 - \)X[, b 42 = y(p 2 2 - pi, b 43 = |m!V(p 2 - 1), c 41 = |(p 2 - l)A, 


M / 


C 42 — - (p2 l)4 > C 43 — 0> Al — 0i fl 51 — G]Q + 
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Pi ~ 1 
Pi 


„ M 
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' p2 1 V, + —MN 
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a 52 2 G]Q, fl 53 


fp? - lA _ M 

i p 2 i r* fl54 2 


1 _ M 


^^L_i]x 2 , £> 53 = \mn 

p2 


f p\ -1 

p 2 


' 2 _ 1 
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fe =i 
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p 3 


* C -=2 


Pi ~ 1 j c = M 
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Pi J 2 


' p2 "I 

Pi 


^r 1 " 


C5 3 — 0, d 32 — 0 


Provided that detailed shapes of functions A, B, C, D, E, H, I , a, and /? are as follows: 
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